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We compute the QCD corrections to the gluon fusion subprocess gg — > 77, which forms an 
important component of the background to the search for a light Higgs boson at the LHC. We 
study the dependence of the improved pp — > 77V background calculation on the factorization and 
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I. INTRODUCTION 

The nature of electroweak symmetry breaking remains 
a mystery, despite decades of theoretical and experimen- 
tal study. In the Standard Model, the masses for the W 
and Z bosons, quarks and charged leptons are all gen- 
erated by the Higgs mechanism. This mechanism leaves 
as its residue the Higgs boson, the one undetected ele- 
mentary particle of the Standard Model, and the only 
scalar Jy . Its properties are completely specified once its 
mass is determined. Alternatives to, or extensions of, the 
Standard Model electroweak symmetry breaking mecha- 
nism typically also include one or more Higgs particles. 
Experiments over the next decade at the Fermilab Teva- 
tron and the CERN Large Hadron Collider (LHC) should 
shed considerable light on electroweak symmetry break- 
ing, in particular by searching for these Higgs bosons, 
and measuring their masses, production cross sections, 
and branching ratios. 

There are good reasons to believe that at least one 
Higgs particle will be fairly light. The Standard Model 
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Higgs boson mass is bounded from above by precision 
electroweak measurements, m# < 196-230 GeV at 95% 
CL Q . In the Minimal Supersymmetric Standard Model 
(MSSM), the lightest Higgs boson is predicted to have 
a mass below about 135 GeV ||; over much of the pa- 
rameter space it has properties reasonably similar to the 
Standard Model Higgs boson. There are also hints of a 
signal in the direct search in e + e~ — > HZ at LEP2, just 
beyond the lower mass limit of 114.1 GeV {§. The corre- 
sponding lower limit on the lightest scalar in the MSSM 
is only 91.0 GeV Q, because the HZZ coupling can be 
suppressed in some regions of parameter space. 

Run II of the Tevatron can exclude Standard Model 
Higgs masses up to « 180 GeV with 15 fb _1 per experi- 
ment. However, at this integrated luminosity a 5<7 discov- 
ery will be difficult to obtain || for a mass much beyond 
the LEP2 limit. Also, the Higgs decay modes relevant for 
searches at the Tevatron, H — ► bb and W + W~, do not 
lend themselves to a precise measurement of the Higgs 
mass. The LHC will completely cover the low mass region 
preferred by precision electroweak fits and the MSSM, as 
well as much higher masses. For run < 140 GeV, the 
most important mode involves production via gluon fu- 
sion, gg — > H, followed by the rare decay into two pho- 
tons, H — > 77 ft 0|. Although this mode has a very 
large continuum 77 background |1 , the narrow width of 
the Higgs boson, combined with the mass resolution of 
order 1% achievable in the LHC detectors, allows one to 
measure the background experimentally and subtract it 
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from a putative signal peak |T^, O, [T^, [l|| . 

For the same mass range of m,H < 140 GeV at the 
LHC, Higgs production via weak boson fusion, qq' — > 
qq'W + W~ — > qq'H, followed by (virtual) weak boson 
decay, H — ► W + W~ — > e ± /i T ^T, is also promising, even 
for Higgs masses as low as the LEP2 limit H|. On the 
other hand, a mass determination from this mode, or 
from weak boson fusion followed by H — > t + t~ [[L5[ , 
cannot compete with the 77 mode, although these modes 
certainly offer very useful branching ratio information. 

The purpose of this paper is to provide an improved 
calculation of the prompt (i.e. not from hadron decay) 
di-photon background to Higgs production at the LHC, 
in particular by computing QCD corrections to the gluon 
fusion subprocess. Although the background will be mea- 
sured at the LHC, it is still useful to have a robust the- 
oretical prediction in order to help validate the quanti- 
tative understanding of detector performance. Perhaps 
more importantly, the theoretical prediction can be used 
to systematically study the dependence of the signal rel- 
ative to the background on various kinematic cuts, pro- 
viding information which can be used to optimize Higgs 
search strategies. As a side benefit, one can improve the 
predictions for a variety of di-photon distributions that 
can be measured. 

The process pp — > 77A proceeds at lowest order via 
the quark annihilation subprocess qq — ► 77, which is 
independent of the strong coupling a s . The next-to- 
lcading-order (NLO) corrections to this subprocess have 
been incorporated into a number of Monte Carlo pro- 
grams Jig, |l7|) . However, the gluon distribution in the 
proton becomes very large at small x, making formally 
higher order corrections involving gluon initial states very 
significant for the production of low-mass systems at the 
LHC. Gluons can fuse to photon pairs through one-loop 
quark box diagrams such as the one shown in fig. |l](a). 
The order a s contribution to pp — ► 77 A from gg — ► 77 
is indeed comparable to the leading-order quark annihi- 
lation contribution |[ |l6| |l|, ]T||. (In fact, the NLO 
correction to quark annihilation, including the qg — » 775 
amplitude, can also be as large as either of these terms.) 
Hence, to reduce the uncertainty on the total 77 pro- 
duction rate, we have computed the contributions of the 
99 ~~ * 77 subprocess at its next-to-leading-order, which 
we shall call "NLO" , even though it is formally N 3 LO as 
far as the whole process pp — > 77 A is concerned. 

The NLO gluon fusion computation has two matrix- 
element ingredients: 

• The virtual corrections to gg — > 77, involving two- 
loop diagrams such as the one in fig. |l](b), which 
were computed recently ]19[| u sing the integration 
methods developed in refs. ]20[ . 

• The effects of gluon bremsstrahlung, through the 
one-loop amplitude for gg — > 773, including pen- 
tagon diagrams as depicted in fig. |l](c). The gg — > 
773 amplitude can be obtained from the one-loop 
five-gluon matrix elements pl[| by summing over 



permutations of the external legs |22|, |23j . 

Both the virtual and real corrections have been evaluated 
in the limit of vanishing quark masses. In the range of di- 
photon invariant masses relevant for the Standard Model 
and MSSM Higgs searches, 90-150 GeV, this is an excel- 
lent approximation. The masses of the five light quarks 
are all much less than the scale of the process. The top 
quark contribution is negligible until the invariant mass 
approaches 2m t ~ 350 GeV; at 150 GeV, it is still well 
under one percent of the total quark loop contribution. 
The virtual and real corrections are separately infrared 
divergent. We have used the dipole formalism p4j to 
combine them into a finite result, in a numerical program 
that can compute general kinematic distributions. 

Prompt photons are not only produced directly in hard 
processes, but also via fragmentation from quarks and 
gluons. As discussed in ref. |l^| , even though the separa- 
tion into direct and fragmentation contributions is some- 
what arbitrary, it is still very useful to track the pieces 
separately. The fragmentation processes occur at low j?t 
with respect to a neighboring jet. They can only be com- 
puted with the aid of nonperturbative information, in 
the form of the quark and gluon fragmentation functions 
to a photon, -D 7 /( 9 or g)i z i^F)- Here z is the fractional 
collinear momentum carried by the photon, measured rel- 
ative to the momentum of the parton which fragments 
into it, and /j,f is the factorization scale used to separate 
the hard and soft processes. The inclusive di-photon pro- 
duction rate at the LHC is actually dominated by the sin- 
gle fragmentation process, e.g. the partonic subprocess 
qg — > 7(7, followed by the fragmentation q — > 7A. Even 
the double fragmentation process, e.g. gg — > qq followed 
by the fragmentations q — > 7 A" and q — > 7A, can ex- 
ceed the direct contribution in the 80-140 GeV range for 
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However, fragmentation contributions can be effi- 
ciently suppressed by photon isolation cuts p7| . Such 
cuts are mandatory in order to suppress the very large 
reducible experimental background where photons are 
faked by jets, or more generally by hadrons. In par- 
ticular, 7r°s at large pt decay into two nearly collinear 
photons, which can be difficult to distinguish from a sin- 
gle photon. The standard method for defining an isolated 
photon is to first draw a circle of radius R in the plane of 
pseudorapidity, r) = lntan(0/2), and azimuthal angle, <fi, 
centered on the photon candidate. The amount of trans- 
verse hadronic energy in this circle, or cone, is required 
to be less than some specified amount, fixmax- Here, R 
and Et max may in principle be varied independently. Al- 
though photon isolation is improved by increasing R and 
decreasing i?Tmax, this cannot be done indefinitely, for 
both theoretical and experimental reasons. Theoretically, 
it is not infrared-safe to forbid all gluon radiation in a fi- 
nite patch of phase space, so any prediction would have 
large uncontrolled corrections. Experimentally, fluctua- 
tions in the number of soft hadrons from the hard scat- 
tering, the underlying event, and other minimum bias 
events in the same bunch crossing, plus detector noise, 
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FIG. 1: Sample quark loop diagrams contributing to pp — > 77X, which are computed in this paper: (a) the leading order 
gluon fusion subprocess gg — > 77, (b) the virtual correction to this subprocess, and (c) the radiative process gg — > 773. 



impose a lower limit on the £/rmax that can be required 
for isolation, at a given R. 

A typical choice of isolation criteria in past LHC stud- 
ies has been R = 0.4 and £?Tmax = 5 or 15 GeV 0. 
The experimental optimization of these variables has 
often been made with the suppression of the huge re- 
ducible background from jet fakes as a primary crite- 
rion STL |l^, [l| . This criterion is very understandable 
in the light of how poorly this background is understood; 
it depends on the tails of distributions, such as the very 
hard (z w 1) tail of parton fragmentation to tt°s p5| . 
Nevertheless, it is estimated that this background can 
be reduced to the order of 10-20% of the irreducible 77 
background, so that one should try to optimize with re- 
spect to the latter background as well. We shall inves- 
tigate the behavior of the irreducible 77 background, as 
well as the Higgs signal, as the cone size R is increased, 
while also increasing the transverse energy allowed into 
the cone. 

An alternative "smooth" cone isolation criterion has 
been proposed by Frixione pq| , in which a continuous 
set of cones with r < R are defined, and the transverse 
energy permitted inside r, E^(r), decreases to zero as r 
does. This criterion has the theoretical advantage that it 
entirely suppresses the more poorly known fragmentation 
contribution, while still being infrared safe. Experimen- 
tally, however, this theoretical ideal may be difficult to 
achieve with a detector of finite granularity, and taking 
into account the transverse extent of the photon's elec- 
tromagnetic shower [^3|. Nevertheless, we shall study 
the behavior of Higgs signal and background for a few 
versions of the smooth cone criterion as well. 

The fragmentation contributions are technically in- 
volved to compute at NLO. Their implementation re- 
quires, for example, all the one-loop four-parton matrix 
elements and tree-level five-parton matrix elements, as 
well as convolution with fragmentation functions. Fortu- 
nately, a flexible program is available, DIPHOX (l7], p5[ , 
which incorporates at NLO the quark annihilation direct 
subprocess, plus single and double fragmentation. It also 
contains the gluon fusion subprocess at leading order. 
We have used DIPHOX to produce all the non-gluon fusion 
contributions to pp — > 77X for the case of standard cone 
isolation. These were then combined with the results of 



our NLO implementation of gluon fusion. (We have also 
checked the DIPHOX implementation of the leading order 
gluon fusion contribution against ours, and they agree 
perfectly.) For the case of smooth cone isolation, where 
fragmentation contributions are absent, we have imple- 
mented the qq — > 77 direct subprocess at NLO ourselves, 
as well as the gluon fusion subprocess. 

Since we are computing only a selection of the contri- 
butions to pp — > "/"/X at order a 2 and , we should make 
a few remarks about contributions we have omitted, at 
least those appearing at order a 2 . Figure || indicates a 
few such contributions. Perhaps most worrisome at first 
glance is the tree-level cross section for gg — > "/"/qq, as 
shown in fig. ||(a). This subprocess can utilize the large 
gluon-gluon luminosity, is only of order a 2 s , and has not 
yet been computed to our knowledge. In addition, it can 
be enhanced by 7(7 collinear singularities if the isolation 
cut is not too severe. It has a double collinear singu- 
larity related to the g — > qq splitting process, so that 
this contribution should depend on the factorization scale 
through terms proportional to a 2 s (\n. lip)™ with n = 1,2. 
Although this contribution alone does not cancel the full 
dependence on fxp to order a 2 , it does cancel the terms 
in da / d(\a. Hf) that are enhanced by two factors of the 
gluon density and are potentially the largest. Therefore, 
it is reasonable to consider the gg — > 77^ contributions 
as a first approximation to the full NNLO calculation. A 
study of the similar subprocesses, gg — > WVqq' where 
V = 7 or Z, was recently carried out, with the conclu- 
sion that their effects are actually quite small, 5% or less 
(and negative) , with respect to the lower-order contribu- 
tions of the qq and qg initial states pTj. Assuming that 
these results also hold for the 77 final state, our neglect 
of gg — > "/"/qq should be tolerable. 

Other order a 2 corrections to the quark annihilation 
subprocess include the two-loop virtual corrections shown 
in fig. |^(b) , which also have been computed recently! 28) , 
and real corrections, such as from qg — > "/"/gq in fig. E(c). 
The numerical implementation of these corrections in- 
volve doubly unresolved parton configurations and re- 
main to be performed. Because the qg — > 775 amplitude 
provides such a large NLO correction to the quark an- 
nihilation subprocess, the order a 2 contributions arising 
from the qg initial state could be reasonably important. 



4 




FIG. 2: Sample diagrams for contributions to pp — > 77X which are not treated in this paper, although they are of order a? a : 
(a) the tree-level subprocess gg — > jjqq, (b) doubly virtual correction to qq — > 77, (c) doubly real correction to qq — > 77, of 
the form — > ^yygq, and (d) the process qg — > 775 wa a quark box. 



The qg initial-state contributions may be more tractable 
than the full order a 2 corrections to the quark annihila- 
tion subprocess, because of fewer soft gluon singularities. 

There are actually three other ways for the quark box 
responsible for gg — > 77 to contribute to the pp — > 77 A 
cross section at order a 2 . The quark box can be incor- 
porated into the one-loop qg — > 77(7 amplitude shown 
in fig. §(d), which can then interfere with the tree- level 
qg — ► 775 amplitude. However, such a contribution only 
gets to utilize one gluon distribution function at small x, 
not two. It is likely to be smaller than the other order a 2 
qg initial-state contributions just mentioned, given that 
it lacks both initial- and final-state singularities. (The 
loop and tree amplitudes each have square-root collincar 
singularities, but they appear in different collincar lim- 
its, so that the phase space integral of the amplitude 
product remains finite.) The quark box contribution to 
the crossed process qq — > 773 is also infrared finite, and 
lacks any small x gluon enhancement, so it should be 
even smaller than the qg — ► 775 case. Finally, there is 
a two-loop virtual correction to qq — > 77 containing the 
quark box (not shown here) , which can interfere with the 
tree amplitude at order a 2 . This particular correction is 
also infrared finite, and does not benefit from any small 
x gluon distribution. Numerical evaluation of the expres- 
sion in ref. shows that its magnitude never exceeds 
0.3% of the Born quark annihilation cross section (its sign 
is negative for all relevant, central scattering angles), so 
we are justified in neglecting it. 

In this paper we also investigate the effects of isola- 
tion cuts and other kinematic features of the Higgs sig- 
nal vs. the QCD background. Intuitively, one might ex- 
pect the irreducible background to diminish more rapidly 
than the signal as one imposes more severe photon iso- 
lation requirements for the following reason: The largest 
component of the irreducible background, for typical iso- 
lation cuts, comes from the subprocess qg — > 779, even 
after we include the NLO corrections to gluon fusion. 
The (/(/-initiated subprocess has final-state singularities 
when the quark is collinear with either of the photons, 
which of course require some kind of photon isolation to 
make finite. Thus, the contribution of this subprocess 
can be significantly reduced with tighter isolation cuts. 
On the other hand, the dominant Higgs production pro- 



cess, gg — ► HX, even at higher orders where additional 
partons are radiated, clearly does not give rise to partons 
that are preferentially near the photons. 

There are at least two different ways to make the pho- 
ton isolation more severe. One way is to strengthen the 
isolation cone, for example by increasing its radius R. 
A second approach is to impose a jet veto |l2j in the 
neighborhood of the photon candidate, on top of a "stan- 
dard" isolation requirement. For example, one may for- 
bid events having a jet with > £/rjet whose axis is 
within a distance i?j e t in 4>) space from either pho- 
ton candidate, where i?j G t > R- The two approaches 
are depicted in fig. |3| The small inner cone corresponds 
to a standard photon isolation cone, e.g. with radius 
R = 0.4. The outer cone could correspond to a larger 
isolation cone, e.g. with radius R = 1 or R = 2. Alter- 
natively, it could represent the radius i?j C t within which 
jets, such as the one depicted, are vetoed against. 

At the level of a NLO parton calculation, the jet veto 
approach is not too different from the large cone ap- 
proach, if one chooses i?j e t and -Bxjet to be equal to the 
large cone values of R and i?Tmax, respectively. The rea- 
son is that for the direct contributions the hadronic en- 
ergy is deposited as a single parton. However, the jet veto 
approach is probably preferable for experimental reasons. 
The hadronic energy represented by the final state quark 
in the background process qg — > 775 should typically be 
deposited in the form of a jet. In contrast, an equivalent 
amount of energy originating from either soft hadrons 
from overlapping events, or detector noise, will usually be 
uniformly distributed in (??,</>), and so will be less likely 
to form a jet. Thus when a signal event is accompanied 
by such energy, it is less likely to be discarded in the case 
of a jet veto, compared with the case of a standard or 
smooth isolation cone of large size. 

Another feature distinguishing the two photons from 
Higgs decay from those produced directly is their angu- 
lar distribution. The qq- and (/(/-initiated components of 
the background have i-channel fermion exchanges that 
tend to yield photons toward smaller scattering angles in 
the 77 center-of-mass (for a fixed M 77 ), compared with 
the decay of the Higgs boson, which is isotropic in cos 9* . 
Hence, we will study signal and background distributions 
in y* = (2/(71) — j/(72))/2, where the y{"fi) are the rapidi- 
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ton isolation criteria, either by increasing the size of the 
isolation cone, or by imposing a veto on nearby jets. We 
also compute the distributions in y* and y 11 , and investi- 
gate how they may be used to help distinguish the Higgs 
signal from the background. In section ^ we present our 
conclusions, and the outlook for further improvement of 
our understanding of the 77 background. 




II. OUTLINE OF THE COMPUTATION 



FIG. 3: Illustration of two possible ways to strengthen pho- 
ton isolation, beyond the "standard" isolation represented by 
the inner cone. One can either increase the cone radius to the 
large outer one shown, or one can veto on jets within such a 
radius. 



ties of the photons. This variable serves as a proxy for 
cos 8*, because cos 8* = tanhy* at lowest order. 

The remaining variable which describes the Higgs bo- 
son kinematics at LO is the rapidity of the Higgs boson 
itself. This leads us to consider the signal and back- 
ground distributions of y 77 , the rapidity of the total di- 
photon system. To a large extent this distribution is de- 
termined by the parton luminosities which are involved in 
the production process. Since Higgs production is domi- 
nated by gluon-gluon fusion, whereas the direct 77 back- 
ground gets comparable contributions from all three ini- 
tial states, qq, qg, and gg, there may be a useful difference 
in the y 77 distribution. 

To carry out these studies of the signal as well, we have 
implemented the gluon fusion production of the Standard 
Model Higgs boson at NLO, followed by decay to 77. We 
work in the heavy top quark limit, for which an effective 
Hgg vertex [pj suffices to describe the production pro- 
cess at the low Higgs transverse momenta we consider. 
(We do, however, include the exact mg /mt dependence 
of the LO term. This has been shown [g9| to be an excel- 
lent approximation to the exact NLO cross section |3(]] 
for the Higgs masses with which we are concerned.) The 
total Higgs production cross section has recently become 
available in this limit at NNLO The corrections 

from NLO to NNLO are modest. Here we wish to study 
distributions with kincmatical cuts, which have not yet 
been computed at NNLO. For cuts that select events 
with nonzero transverse momentum of the Higgs boson 
(which we do not explicitly consider in this paper), our 
calculation is only a LO computation; distributions for 
such quantities at NLO are available elsewhere (32|. For 
the branching ratio Br(_ff — > 77) we use the program 
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This paper is organized as follows. In section || we out- 
line the NLO gluon fusion computation. In section |l| wc 
show the effects of the NLO corrections to gg — > 77 on 
the total pp — ► 77X cross section, using kinematic cuts 
appropriate to the Higgs search, and concentrating on 
the comparison with LO and on the scale dependence. 
In section IV we study the effects of varying the pho- 



The gluon fusion subprocess begins at one loop. The 
next-to-leading order corrections to it include two-loop 
four-point amplitudes and one-loop five-point ampli- 
tudes. The calculation of these matrix elements requires 
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a fair amount of work |1S|, [21], |2S 
the amplitudes are available, the rest of the NLO gluon 
fusion computation is in all essential respects identical 
to a conventional NLO computation involving tree-level 
and one-loop amplitudes. For example, the two-loop am- 
plitudes for gg — > 77 have virtual infrared divergences 
with the same structure as those of a typical one-loop 
amplitude. Similarly, the vanishing of the gg — > 77 tree 
amplitude causes the singular behavior of the gg — > 773 
one-loop amplitude, as the outgoing gluon becomes soft 
or collinear with an incoming gluon, to have the same 
form as that of a typical tree amplitude with one addi- 
tional radiated gluon. 

The one-loop gluon fusion helicity amplitudes for the 
process 

3(-Pi,-Ai)+3(-p 2) -A 2 ) 7(P3,A 3 )+7(p 4 ,A 4 ), (1) 

in an "all-outgoing" labelling convention for the momenta 
Pi and helicities Ai, are given by 
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where a is the QED coupling at zero momentum transfer, 
a = 1/137.036 . . . and a s (/ifl) is the running QCD cou- 
pling in MS scheme, evaluated at renormalization scale 
fiR. Also, in this formula, are the adjoint gluon color 
indices, Qj are the quark charges in units of e, the ap- 
propriate number of light flavors is Nf = 5, and we have 
suppressed overall phases in the amplitudes. The quan- 
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where s = (pi + p 2 ) 2 , * = {P2 +P3) 2 , " = (pi +P3) 2 , and 
the remaining helicity amplitudes are obtained by parity. 

The dimensionally-regularized and renormalized two- 
loopQCD corrections to these amplitudes can be written 
as B 
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where the number of dimensions is £> = 4 — 2e, and all 
the poles in e are contained in 
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and the number of colors is N = 3. The finite expres- 
sions xxx ( s ' an d -?a L a a a ( s i are presented in 
ref. . ' 3 4 

The radiative process, 

g(-pi, -Ai) + ,9(-p 2 , -A 2 ) 

-* 7(^3, A 3 ) +7(P4,A 4 ) +5(p 5 ,A 5 ), (7) 

begins at order af. The squared matrix element, aver- 
aged over initial hclicities and colors, and summed over 
final ones, and with a 1/2 for identical final-state pho- 
tons, is given by p2, E3| 
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where cr; label the hclicities and momenta of the gluons 
and photons, and COP^ 125 ^ denotes the subset of 12 per- 
mutations of (1,2,3,4,5) that leave 5 fixed and preserve 
the cyclic ordering of (1,2,5), i.e. 1 appears before 2. 



The partial amplitudes ^j^ 2 ' (<r) are those for five-gluon 
scattering via a quark loop given in ref. pl| , but with an 
overall factor of (4ir)~ 2 removed. 

The permutation sum in eq. (^) cancels out all of the 
virtual divergences of the partial amplitudes, and most of 
their singularities as momenta become soft and collincar. 
The only remaining singularities are when the final gluon 
momentum ps becomes soft, or becomes collinear with 
either initial gluon momentum, p\ or pi. In the region, 
for example, where P5 is collinear with p\, with p$ — > 
— (1 — x)pi, the squared matrix element has the limiting 
behavior, 
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and the gg — > 77 process has the kinematics (— xpi) + 
(— P2) - * P3 + P4. The second term involves the spinor 
products (i j) |34| entering the five-point partial ampli- 
tudes. It accounts for nontrivial phase behavior of the 
amplitudes as p§ rotates azimuthally about the p\ di- 
rection. The primed sum is defined analogously to the 
unprimed sum, except that in each complex conjugated 
helicity amplitude the helicity of the gluon with momen- 



M 



(i) 



M (1) * 

1 ( — Ai)A 2 A 3 A4 



turn (— xp\) is flipped, so that 
appears, and the spinor product overall phases in the 
amplitudes should also be included Q. 

Integration over these singular phase space regions can 
be handled by the dipole formalism pi|, with just two 
dipole subtractions, one for each initial gluon. The dipole 
subtraction for gluon 1 is given by eq. (^0|), where x = 
(Pi + P2 + Pa) 2 /{pi + P2) 2 , and the four-point matrix- 
elements are evaluated for boosted kinematics: gluon 1 
is assigned momentum (— xpi), gluon 2 momentum — P2 
still, and the photon momenta pj, j = 3, 4, are set equal 
to 



P»=P»- 2p l} K t K \ K + KT 



where 



(K + Kf 

= 
= 



2pj ■ K - 



K^, (12) 



-Pi 

-XPl 



a a 

P2-P51 
P2 ■ 



(13) 



Subtracting these dipole terms from the gg — > 7717 
squared matrix element, eq. (@), removes the soft and 
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collinear singularities, so that the resulting expression 
can be integrated directly in four dimensions. Further- 
more, the dipole subtraction terms themselves can be in- 
tegrated analytically in D = 4 — 2e dimensions [H; their 
poles in e, combined with the collinear counterterm in 
the MS factorization scheme, cancel against the virtual 
divergence in the interference of M 2 ~}%°v in eq. (|J) with 
Aigg^°v. Thus, by adding and subtracting the dipole 
terms, one obtains an expression for the NLO cross sec- 
tion which is explicitly finite in four dimensions. 

The remaining finite NLO differential cross section for 
PP ~ * 99 ~ * 11 X consists of three terms, which can be 
put in the form 



(1(7 



NLO 



dxi dx 2 dz ,...„/ o / \ 

. g(xi)g(x 2 )6 (z - m^/s) 



x 6(1 - z)da B + da G + da K , 



(14) 



where the three terms are functions of the incoming par- 
ton momenta x\p a and X2Pb and the variable z, and we 
have used s = (x\p a + X2Pb) 2 ■ The term with leading- 
order kinematics ([!]), including the LO answer, can be 
written 



da = do 



LO 



1 



oc s (hr) 



— 2 Re 

2s 



Ea / |2— loop, fin. A/fl — loop* 
J 99^11 J 99^11 



hel. , color 



dT 2 , (15) 



where 



where dT$ is the three-particle Lorentz-invariant phase 
space. Thus it involves the full five-point kinematics. 

We have implemented these three terms in a numeri- 
cal program that allows for different kinematic cuts and 
photon isolation criteria to be applied. The numerical 
integrals are performed by adaptive Monte Carlo sam- 
pling using VEGAS pEfl . For the first and second terms 
the Monte Carlo routine is used to generate the photon 
four-momenta with four-point kinematics, along with the 
appropriate longitudinal boosts of eq. (p7[). For the third 
term the events are generated with five-point kinematics; 
however, treatment of the final-state particles differs for 
the gg — » 77(7 squared matrix element eq. (^|) and the 
two dipole subtractions. Whereas the gg — > 7717 squared 
matrix element is treated as a true three-parton final- 
state, the dipole subtraction for gluon 1 is treated as a 
two-parton final-state with the photon momenta given 
by eq. (12), and similarly for the dipole subtraction for 
gluon 2. Thus, each call to VEGAS in the third term pro- 
duces three distinct kinematic configurations with three 
different weights. The infrared safety of the kinematic 
and isolation cuts ensures the appropriate cancellation 
between the pieces as the gluon becomes soft or collinear. 



III. RESULTS FOR DI-PHOTON 
BACKGROUND 

In this section we present the modifications to the 
PP ~ * llX cross section due to the inclusion of the NLO 
contributions to the gluon fusion subprocess. 



da 



LO 



da {xip a ,x 2 Pb) 



2s 



E \M£Z%\ a dT 2 



hcl. .color 



dT 2 is the two-particle Lorentz-invariant phase space, and 
■M?ji?2y' fln ' refers to just the terms containing F L and 
F SL in eq. (||). The second term also has leading-order 
kinematics, but boosted along the beam axis, with 

da c = ^lzK^(z) 
2tt v ; 

da LO (zx 1 p a ,X2Pb) + da LO (x 1 p a ,zx 2 Pb)] , (17) 
ln((l - z)M 11 /p, F ) 



where 

K 9 ' 9 (z) = 2N 



1-z 



In ; 



1 



1 + *(!-*) 



(1 - z) 2 M 2 
In I iL^! XL 

z p A F 



.(18) 



The third term contains the gg — > 7717 squared matrix 
element (||) , minus the two dipole subtractions mentioned 
above, 



da* = - 
2s 



\M\ 2 ad -\M\ 2 ipolcA -\M\ 2 ipolCi2 



dT 3 , (19) 



A. General remarks 

We impose the following kinematic cuts on the two 
photons, 

pr(7i) > 40 GeV, p T ( l2 ) > 25 GeV, \y{n, 2 )\ < 2.5, 

(20) 

which are essentially those used by the ATLAS and CMS 
detectors in their Higgs search studies |l(], 0. In addi- 
tion, we require each photon to be isolated from hadronic 
energy, according to one of two criteria: 

1. standard cone isolation — the amount of trans- 
verse hadronic energy Et in a cone of radius R — 
\J (Arj) 2 + (A</>) 2 must be less than St max- 

2. smooth cone isolation |26[ — the amount of trans- 
verse hadronic energy Et in all cones of radius r 
with r < R must be less than 



E 



T max 



(r) =7J T (7) e 



1 — cos r 
1 — cos R 



(21) 



Here we shall choose n = 1. 



Note that in an NLO calculation, "hadronic energy" in 
the neighborhood of the photon always amounts to just 
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a single parton, a rather crude approximation to the true 
hadronic background. 

Unless otherwise specified, we evaluate the NLO cross 
sections with the MRST99 set 2 parton distributions |36| , 
with the corresponding value of a s (mz) = 0.1175. For 
comparison purposes, we also present cross sections for 
the leading-order gluon fusion subprocess, convoluted 
with the NLO parton distributions MRST99 set 2. We 
use NLO instead of LO parton distributions here, because 
that a ppr oach was taken in the most complete previous 



study 1 17 , and we wish to highlight differences with re- 
spect to that work. The use of NLO parton distributions 
can also be justified by considering the "LO" gluon fu- 
sion subprocess as a NNLO correction to the entire 77 
production process. 

Our default choices for the renormalization and factor- 
ization scales are 



A*if = M-F = 0.5M, 



77' 



(22) 



as in rcf. JL7|; we also investigate the dependence of the 
results on (1r and [ip. For the fragmentation contribu- 
tions jl?]] , we rely on the NLO set I photon fragmentation 
functions from rcf. |S7|. 



B. Effects of NLO gluon fusion on the 77 
background 

Figure ||(a) shows the contribution of just the gluon fu- 
sion subprocess to pp — > 77 V at the LHC, at its leading 
and next-to-leading orders, for the standard cone photon 
isolation criterion with R = 0.4, Ermax = 15 GeV, and 
for several choices of parton distributions. To help with 
comparisons to the results of ref. 0, we use MRST99 
set 2 as our "default" choice. This set has a somewhat 
larger gluon distribution at large x than MRST99 set 1, 
but the differences with this set, or with CTEQ5M1 @, 
at the smaller x ranges probed here are small compared 
to the NLO corrections, or to the renormalization and 
factorization scale dependence, as we shall see. We also 
plot the LO cross section with the LO CTEQ5L distri- 
butions (using a LO a s (mz) = 0.127) in order to 
compare our NLO cross section with a "true" LO calcu- 
lation. Recently, the more precise HERA data has been 
incorporated into two updated standard sets of distri- 
butions, MRST2001 [|9) and CTEQ6M @. However, 
neither set has a sizable change from its predecessor in 
the quark and gluon distributions for x in the relevant 
range 0.01—0.1 at Q 2 = 10 4 GeV 2 . The new MRST2001 
distribution uses a slightly larger a s (mz) — 0.119, which 
may increase the importance of the gluon fusion subpro- 
cess relative to the qq subprocess by a few percent, but 
overall the effect on the 77 background should be fairly 
small. 

In the absence of an NLO calculation, some experi- 
mental studies had used the K factor (ratio of NLO over 
LO cross section) for Higgs production by gluon fusion 



TABLE I: NLO QCD K factors for 77 Higgs signal and 
gluon fusion background. Both LO and NLO cross sections 
are computed using NLO parton distributions. 



M 77 (GeV) 


^Higgs 




98 


2.92 


1.82 


118 


2.54 


1.61 


138 


2.39 


1.55 



as an estimate of the K factor for gg — > 77. For exam- 
ple, -ft^HWs = was uset ^ f° r a 100 GeV mass Higgs 
boson pj. The reasoning is that both gg — > H and 
gg — > 77 involve production of a colorless system from a 
gluon-gluon initial state. One difference between the two 
processes, however, is that the Hgg coupling receives a 
fairly large short-distance renormalization p0| , from the 
top mass scale, 



Higgs 



1 



2 7T 



1.2, 



(23) 



which has no counterpart in the gg — ► 77 correction. An- 
other difference stems from the different loop momentum 
scales appearing in the real emission diagrams. For the 
gg — > Hg case, the momentum in the loop is dominated 
by the heavy top quark mass, which is taken to be infinite 
in our calculation, while for the gg — > 7717 process the 
quarks in the loop are taken to be massless, so that the 
dominant momentum in the loop is that of the photons 
and the emitted gluon. As a result, the cross section for 
gg — > 7717 falls off more quickly with the emitted gluon 
transverse momentum than that for gg — ► Hg, resulting 
in a smaller real-emission contribution to the total NLO 
cross section. 

In comparing signal and background K factors, it is 
of course useful to impose the same set of cuts on the 
photons in each case. In table | we list K factors for 
both the pp — > gg — > HX — > 77V Higgs production 
cross section and the pp — > gg — ► 77V gluon fusion back- 
ground, for three representative choices of Higgs mass. 
(To be precise, the Higgs K factor includes the subpro- 
cesses qg — -> Hq and qq — > Hg at NLO; removing them 
decreases -R'mggs by roughly 5% at run = 118 GeV.) We 
take /ifj ; = jip — 0.5M 77 and impose the same photon ac- 
ceptance and isolation cuts as in fig. ^(a), for both signal 
and background. (With both sets of cuts removed, each 
K factor is about 10% smaller at M 77 = 118 GeV, but 
their ratio is stable to a few percent.) We define the "LO 
cross section" entering the K factor using NLO, rather 
than LO, parton distributions. This convention results 
in larger K factors than the more standard convention, 
as can easily be seen by comparing the LO cross sec- 
tions using the CTEQ5M1 and CTEQ5L distributions in 
fig. f|(a). In any case, the important point is that the K 
factors for the gluon fusion component of the di-photon 
background are significantly smaller than the K factors 
for the Higgs signal, even after accounting for the short- 
distance contribution (|23|) to the latter. This difference 
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FIG. 4: (a) Contribution of the gluon fusion subprocess to pp — > 77X at the LHC, at leading order (lower four), and at 
NLO (upper three), using various parton distributions, (b) Total pp — > 77X production at NLO, including NLO qq — > 77 
and fragmentation contributions, with the gluon fusion subprocess treated at LO (dashed) and at NLO (solid). MRST99 set 
2 partons are used in (b). Contributions not involving gluon fusion into photons are obtained from DIPHDX. Both sets of plots 
are for fiR = /ip = 0.5M 77 , and a standard photon isolation criterion with R — 0.4, Eimsx = 15 GeV. Statistical errors from 
numerical integration are shown. 




M (GeV) M (GeV) 



FIG. 5: Scale dependence of (a) the gluon fusion subprocess contribution to pp —> 77X, and (b) the total pp — > 77X production 
cross section, for standard photon isolation with R = 0.4, i?Tmax = 15 GeV. In both plots, the bands represent the result of 
varying fin and (if over the square region 0.5M 77 < (ir, fj,p < 2M 77 . The dashed (solid) hatched band corresponds to including 
the gluon fusion subprocess at LO (NLO). For the leading order band in (a) only, the LO CTEQ5L parton distributions were 
used; otherwise the NLO MRST99 set 2 distributions were employed. 



appears to be due to the relatively smaller real-emission 
contribution to the background. 

In fig. ||(b) the effects of computing the gluon fusion 
subprocess at NLO are shown, for the total NLO pp — > 
77X production rate, i.e. including also the qq — > 77 
and fragmentation contributions at NLO obtained from 
DIPHQX 0. As in fig. |(a), fi R = fi F = M 77 /2 and the 
isolation criterion R = 0.4, i?Tmax = 15 GeV is used. 
(The lower histogram in fig. ^(b), where the gluon fusion 
subprocess is treated at LO, corresponds to the result 



in fig. 13 of rcf. [O, except that R = 0.4, Ermax = 
5 GeV is used in that plot.) The increase in the total 
irreducible 77 background which results from replacing 
the LO gluon fusion quark box by the NLO computation 
is a modest one, except at the lowest invariant masses 
relevant only for non-Standard-Model Higgs searches. In 
fig. §(b), the increase ranges from 27% at M 77 = 82 GeV, 
to 10% at 100 GeV, and only 6% at 138 GeV. For the 
most interesting mass range for the Higgs boson in this 
channel, 115 GeV < tuh < 140 GeV, the overall effect 
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FIG. 6: Scale dependence of (a) the gluon fusion subprocess contribution to pp — > 77X, and (b) the total pp — ► 77 V production 
cross section. The cuts are identical to the ones in fig. 0. The parameters are identical to the ones used in the corresponding plots 
in fig. |f| except that the renormalization and factorization scales are set equal and varied between 0.5M 77 < (j,r = < 2M 77 . 



on the square root of the background is under 5%. The 
larger increase at smaller M 77 simply reflects the fact 
that the LO contribution vanishes, due to the kinematic 
cuts, as M 77 — > 80 GeV. This feature is seen most visibly 
in fig. |(a). 

In fig. [|, the dependence of the 77 background on the 
renormalization scale and factorization scale /if is il- 
lustrated by varying them independently over the square 
region 0.5Af 77 < hr,^f < 2M 77 . Figure ||(a) shows 
the variation for the gluon fusion subprocess contribu- 
tion alone, while fig. g|(b) shows the variation for the 
total production rate, treating the qq — > 77 and fragmen- 
tation contributions at NLO. The same photon isolation 
criterion is used as in fig. ||. In fig. ||(a), the leading- 
order (dashed hatched) band is computed using the LO 
parton distribution CTEQ5L, which is a bit more ap- 
propriate when considering this subprocess in isolation. 
In all cases, the maximum cross section in the band 
at a given M 77 comes from setting /j.r — 0.5A/ 77 and 
fip = 2M 77 , while the minimum cross section comes 
from setting fi R ; = 2M 77 and //> = 0.5M 77 . Allowing 
independent variations for [ar and /zf results in NLO 
bands which are not appreciably narrower than the LO 
bands. In contrast, varying \xr and \xf together, i.e. 
fJ'R = fJ-F = X^ti with 0.5 < x < 2, as is more con- 
ventional, leads to much less scale variation for the gluon 
fusion subprocess at NLO than at LO, as shown in fig. ||. 
(These general features are also qualitatively present in 
the Higgs production cross section as well, although the 
larger NLO K factor leads to stronger renormalization 
scale dependence in that case; see e.g. ref. pi].) The 
considerable improvement in the scale variation for the 
gg contribution depicted in fig. ||(a) is diluted in fig. |^(b), 
where the contributions with quark initial states and 
fragmentation are added in. 

In conclusion, the NLO corrections to the gg — > 77 



subprocess have a modest effect on the total irreducible 
di-photon background to the Higgs search. Thus this 
subprocess can be considered to be under adequate the- 
oretical control. 



IV. STATISTICAL SIGNIFICANCE OF HIGGS 
SIGNAL 

In this section we investigate the kinematic features 
of the Higgs signal and background, starting with pho- 
ton isolation criteria. To facilitate this study we con- 
sider a crude approximation to an experimental analysis 
at the LHC. We assume a Higgs mass of 118 GeV, and 
we count the number of events in mass bins of 4 GeV 
for 30 fb _1 of integrated luminosity, corresponding to 3 
years of running at low luminosity, C = 10 33 cm _2 s _1 . 
We note that this choice of mass bin is slightly larger 
than the optimized mass bins of 2.74 and 3.44 GeV used 
in the ATLAS study of ref. We also include an ef- 
ficiency factor of 0.57 for both signal and background, 
corresponding to the combination of 0.81 per 7 for 7/ jet 
identification and 0.87 for fiducial cuts (mainly the tran- 
sition between barrel and endcap) found in that analy- 
sis. Finally, we include a reducible background of 20% of 
the 77 continuum background, which we assume is pos- 
sible after the 7 /jet identification p2[ . For this analysis 
we take the efficiency factors and the percentage of re- 
ducible background as independent of the isolation cuts; 
to investigate this further would require a more serious 
experimental analysis, beyond the scope of this work. 

In computing the statistical significance we ignore in- 
terference between the Higgs signal and the background. 
In the Standard Model, the interference terms are on the 
order of a few percent of the Higgs signal E^] , and do not 
significantly alter any of our conclusions. The small size 
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of the interference is due mostly to the extreme sharpness 
of the underlying Higgs resonance which, before smear- 
ing with the detector resolution, gives a peak in the cross 
section rising about a factor of a hundred over the back- 
ground. Hence the interference contribution should not 
be more than about 20 percent of the signal. However, 
it is less than this because the primary Higgs production 
mode is via gluon fusion, so only the gg — > 77 com- 
ponent of the background can interfere. Also, because 
the experimental width is much greater than the intrin- 
sic width, only the integral in M 77 across the lineshape 
is observable. This integral vanishes unless there is a 
relative phase between the production {gg — ► H), de- 
cay (H — > 77), and background (gg — > 77) amplitudes. 
The phase happens to vanish, up to small quark mass 
effects, when the background amplitude (for identical- 
helicity photons) is evaluated at one loop. The extra 
power of a s in the two-loop amplitude then provides an 
additional suppression factor. 



A. Effects of varying photon isolation criteria 

We first consider the effects on both signal and back- 
ground of varying the photon isolation criteria, before 
turning in section IV B to the phenomenologically more 



viable method of using a jet veto to impose more strin- 
gent cuts. As mentioned in the introduction, photon iso- 
lation can be achieved by either a standard or a smooth 



cone criterion. In section III B we presented cross section 
results for the standard cone criterion with R = 0.4 and 

Sp max = 

15 GeV, values typical to previous analyses. 
Now we shall investigate how the 77 background varies, 
relative to the H — > 77 signal, as we change the isola- 
tion criteria. In particular, we would like to determine 
the parton-level statistical significance of the signal as a 
function of photon isolation. 

Figure |?](a) shows how the pp — > 77 A production rate 
at the LHC depends on the parameters R and Sr m ax of 
the standard cone isolation definition, while Figure f7|(b) 
presents analogous information for the smooth cone cri- 
terion. As the isolation becomes more severe, i.e. R 
is increased or St max or e are decreased, the direct 
pp — > 77 A background becomes more suppressed. The 
large sensitivity to these parameters is indicative of the 
57 collinear singularity in the NLO qq — > 77A cross sec- 
tion. Since the QCD radiation in Higgs production has 
no such singularity, it should have no correlation with the 
photon directions, and therefore it should be less sensi- 
tive to the isolation criterion. 

To see this more clearly, it is instructive to plot indi- 
vidually the various subprocess contributions to the 77 
background. Figure ||(a) and (b), plotted for two dif- 
ferent smooth cone isolation criteria, show how the qg 
component is reduced relative to qq and gg as the iso- 
lation requirement is made more severe. For instance, 
in the bin centered at M 77 = 118 GeV, the qg compo- 
nent decreases by 36%, in going from R — 0.4, e = 1 to 



TABLE II: Number of signal and background events and the 
statistical significance for pp — * HX — > 77X in a bin 116 
GeV < M 77 < 120 GeV for 30 fb" 1 integrated luminosity 
using standard cone isolation. The Higgs mass is taken to be 
rriH = 118 GeV. Other experimental assumptions are given 
in the text. 



(7?,£ T max in GeV) 


S 


B 


s/Vb 


(0.4,15) 


993 


20,400 


7.0 


(0.4,5) 


980 


19,000 


7.1 


(1,30) 


979 


20,600 


6.8 


(1,15) 


952 


16,900 


7.3 


(2,30) 


896 


16,400 


7.0 


(2,15) 


789 


11,000 


7.5 


TABLE III: Number of 


signal and background events and 


the statistical significance 


for pp — > 


HX — > 77V in 


a bin 116 


GeV < M 77 < 120 GeV for 30 ft 


1 integrated luminosity 


using smooth cone isolation. The Higgs mass is taken to be 


inn = 118 GeV. Other experimental assumptions 


are given 


in the text. 








(R,e) S 




B 


S/s/B 


(0.4,2) 993 




22,000 


6.7 


(0.4,1) 992 




20,800 


6.9 


(0.4,0.5) 985 




20,000 


7.0 


(1,2) 969 




18,100 


7.2 


(1,1) 948 




16,700 


7.3 


(1,.5) 915 




15,400 


7.4 


(2,2) 893 




14,700 


7.4 


(2,1) 806 




12,300 


7.3 


(2,.5) 685 




9,800 


6.9 



R = 1, e=l, while the qq and gg components each only 
decrease by about 4%. The smooth cone criterion was 
used to simplify the discussion, since there are no frag- 
mentation contributions; however, the results are quali- 
tatively similar for the standard cone isolation. 

We give the number of Higgs signal (S) and back- 
ground (B) events and the statistical significance 
(S/y/B) in the bin 116 GeV < M 77 < 120 GeV for 
the different choices of standard cone parameters in ta- 
ble |n|, an d for the choices of smooth cone parameters 
in table [II. For the standard cone, we find a statis- 
tical significance of 7.3 at R = 1, £r m ax = 15 GeV, 



for a modest gain of about 4% over the value of 7.0 
at R — 0.4, fixmax = 15 GeV. The significance can 
be increased further to 7.5 by the very severe cut of 
R = 2, St max = 15 GeV, for a gain of 7% over the 
value at R — 0.4, St max = 15 GeV. For the smooth 
cone, the statistical significance appears to have a maxi- 
mum of about 7.4 for R — 1, e — 0.5 and R = 2, e = 2. 
It is clear from these results that for either the smooth 
or standard cones the statistical significance depends on 
the isolation cuts only rather weakly. 
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FIG. 7: Dependence of pp — > 77X at the LHC on photon isolation cuts, for (a) a set of standard cone isolation parameters, 
R and Bimai, and (b) a set of smooth cone isolation parameters, R and e. All plots are for MRST99 set 2 partons, and 
fJ-R = = 0.5M 77 . 
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FIG. 8: Contributions of the different subprocesses to pp — > 77X at the LHC, as a function of di-photon invariant mass, for 
two choices of smooth cone isolation parameters: (a) R = 0.4, e = 1, and (b) R = 1, e = 1. The plots are computed at NLO, 
for MRST99 set 2 partons, with acceptance cuts (Efj), and (ir — /if = 0.5M ly . 



The smallest cone size, R — 0.4, used in tables |l| 
and [II, is the standard cone size used in previous stud- 
ies. Recently it was observed [^3) that logarithms of the 
form a s (/j)\n(l/R 2 ) Q can invalidate an NLO calcula- 
tion of prompt photon production: For R = 0.1 the NLO 
single-photon cross section with isolation was larger than 
the cross section with no isolation, which is clearly an un- 
physical result. It is not yet known precisely how small 
R can be taken before the NLO calculation begins to 
break down. However, for R = 0.4, the ln(l/i? 2 ) factor 
is 2.5 times smaller than it is for the pathological case 
of R = 0.1. Also, the physical scale in the di-photon 
invariant masses relevant for the LHC Higgs search is 
significantly higher than the single-photon pt(j) = 15 
GeV case studied in ref. |H, rendering a s (fi) smaller as 



well. Finally, we note that the large logarithms do not 
arise from the gluon fusion contributions, because there 
is no <?7 collinear singularity. Hence the effect of the log- 
arithms in the single-photon case (where gluon fusion is 
not important and was not included) is diluted somewhat 
in the di-photon case by the gluon fusion contribution. 
Nevertheless, further study of this situation, including 
possibly resummation of the logarithms, could be help- 
ful. 

A more critical issue for this analysis is that the most 
severe isolation parameters may not be phenomenolog- 
ically viable, for both the theoretical and experimental 
reasons mentioned in the introduction. Theoretically, it 
is not infrared safe to forbid all gluon radiation into any 
finite region of phase space. If the isolation criteria ap- 
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proach this limit, the perturbative predictions become 
subject to large corrections and therefore become unre- 
liable H @. After all, two R = 2 cones can cover 
most of the (f?,</>) plane within the detector acceptance, 
and £/rmax = 15 GeV is not a lot of energy at the LHC. 
On the experimental side, the efficiency for collecting sig- 
nal events may decline for reasons that are absent from 
the NLO Monte Carlo. Instrumental (calorimeter) noise, 
pile up, and energy deposition from the underlying event 
plus overlapping minimum bias collisions, all contribute 
an average energy in a cone which scales roughly as the 
area of the cone. Thus one might expect that when R is 
increased, one should also incre£LSG -St max 5 

roughly like 

R 2 . For the R — 0.4 cone typically used, pile up and un- 
derlying events start to saturate the cone at -&rmax ~ 2.5 
GeV [|l3| [HI . For R = 2, saturation would most likely 
be occurring at i?Tmax = 15 GeV, perhaps even at 30 
GeV. Another potential problem is that as one varies the 
cuts to reduce the irreducible contributions, one must be 
sure that the reducible contributions do not get larger, 
undoing the improvement. For example, we have raised 
the total transverse energy allowed near the photon, in 
going from R = 0.4, i?Tmax = 15 GeV to, say, R = 2, 
-^Tinax = 30 GeV, and this may allow more of the re- 
ducible background to enter. This question could be ad- 
dressed by studies along the lines of ref. [ ps| . In any 
case, it is clear that a phenomenologically more sensible 
method for rejecting events with hadronic energy near 
the photons is required. 



B. Jet veto 

As mentioned in the introduction, a veto on nearby 
jets offers another way to suppress the QCD background, 
in particular the qg — > 77(7 process [ pL2| . At the NLO 
parton level, at least for direct processes, it corresponds 
closely to increasing the size of the cone. However, be- 
cause transverse energy is being forbidden into a smaller 
area (the jet cone size), for the same amount of suppres- 
sion at NLO, the jet veto is a more infrared safe criterion, 
and it should also have better experimental properties 
(less loss of signal due to noise, overlapping events, etc.). 

Jet vetoes have been considered previously in search 
strategies for other Higgs decay modes, particularly H — > 
W+W~ e±^ T igggg, In those cases, 
typically a general veto is applied on all jets in the detec- 
tor acceptance with above a certain value. For the 
77 mode, we would only like to veto on jets "close" to 
the photon candidates. Such nearby jets are more likely 
to come from the qg — > 779 subprocess, because of the 
final state 97 collinear singularity, than from the Higgs 
production process gg — > HX. On the other hand, be- 
cause the gluon is in a larger color representation than the 
quark, g<?-initiated production of a color singlet object 
tends to be jettier overall than production initiated by 
qq or qg initial states. In fact, cuts requiring a minimum 
transverse momentum of the 77 pair, Qt > 30 GeV, have 



PP ~* 77X Cross Section at LHC 




M 77 (GeV) 

FIG. 9: Effect on the background process pp — > 77V of veto- 
ing on jets with Et > 15 GeV within Rj at — 2 of either photon 
candidate. The veto is on top of standard cone isolation with 
R = 0.4, £ T max = 5 GeV. MRST99 set 2 partons are used, 
with the default acceptance cuts (Eo|) and scale choice (|22[). 



been proposed to take advantage of this fact |^2|, ^3], Q , 
and enhance the signal. So only jets "sufficiently" near 
a photon candidate should be vetoed. 

We implement the jet veto on top of a standard photon 
isolation cone, represented by the inner cone in fig. ||. We 
require that there is no jet with a transverse energy Et > 
-^Tjot within a radius FLct of the photon, represented by 
the outer cone in fig. [| We do not include hadronic 
energy inside the inner cone in defining this jet. Then 
the results at the NLO parton level do not depend on 
the cone size i??° t nc used in the jet algorithm, but for 
definiteness we suppose i?j™ t nc = 0.7. 

As an example of the jet veto suppression, Figure ^ 
shows the background suppression obtained for a jet veto 
using i?j Ct = 2 and -Sr/jct = 15 GeV on top of a standard 
isolation cone with R = 0.4 and Eimax = 5 GeV. For 
this standard isolation cone the fragmentation contribu- 
tion is rather small, amounting to about 10 percent of 
the total. This simplifies the calculation of the jet veto 
since we can ignore the action of the jet veto on the small 
fragmentation part. For the direct piece at NLO, a jet to 
be vetoed amounts to a lone parton with transverse en- 
ergy Et > -^Tjct = 15 GeV between the inner and outer 
cones 0.4 < R < 2. By ignoring the jet veto rejection 
of the fragmentation term, the background is overesti- 
mated by a few percent. In this approximation, with 
m H = 118 GeV, the bin 116 GeV < M 77 < 120 GeV has 
776 signal events and 12,600 background events, leading 
to a statistical significance of S/\f~B — 6.9. Even though 
the background drops from 19,000 events for the stan- 
dard isolation case with R = 0.4 and i?Tmax = 5 GeV to 
12,600 events when the jet veto is included, the statisti- 
cal significance is essentially unchanged compared to this 
case. This illustrates the rather disappointing insensitiv- 
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where the sum is over the bins in y* 
be compared with a significance of 



0.75 
y* 



FIG. 10: The angular distribution as a function of y* for 
116 GeV < M 77 < 120 GeV. The renormalization and factor- 
ization scales are fiR = fip = 0.5M 77 and the smooth cone 
isolation cuts are R = l,e = 1. The Higgs mass is taken 
to be rriH = 118 GeV. The peaks near y* = 0.94 occur in 
a region where the NLO calculation breaks down and is not 
trustworthy. 



ity of the statistical significance to the presence of the jet 
veto. 



C. Kinematic distributions of signal and 
background photons 

The situation can be improved somewhat by including 
information from the photon angular distribution. Since 
the Higgs boson is a scalar, its decay to two photons is 
isotropic in its rest frame. In contrast, the 77 background 
processes tend to be more peaked toward the beam axis. 
Thus, the angular distribution of the photons can help 
separate the signal from the background. 



Figure 10 shows the normalized distribution in the 
di-photon rapidity difference, y* = (2/(71) — y(72))/2. 
This variable is convenient because it is simple to de- 
termine experimentally, and at lowest order it is related 
to the center-of-mass scattering angle 9* for qq — > 77 or 
99 ~ * 11 by cos#* = tanhj/*. The renormalization and 
factorization scales are set to our default values (|||), and 
only events in the mass bin 116 GeV < M 77 < 120 GeV 
are kept. The smooth cone isolation is used with pa- 
rameters R = 1 and e = 1; similar distributions are ob- 
tained using a standard cone isolation. As can be seen in 
fig. [l(] the angular distribution of the Higgs signal events 
is rather different from the background. We can esti- 
mate the significance that could be obtained by using a 
maximum likelihood function with this information to be 




(E^) 1/2 



7.3 



This number is to 



(25) 



= 7.7 



(24) 



without using the angular information. The 4% relative 
improvement in significance should also hold roughly for 
y* distributions constructed using a standard cone isola- 
tion. 

An interesting feature in fig. |l^ are the peaks in both 
the signal and background in the bins near 0.90-1.00. 
These peaks are attributable to a breakdown of the NLO 
approximation near the LO kinematic boundary in y*, 
whose location is dictated by the pt(ii) > 40 GeV 
cut @ and A/ 77 w 118 GeV. At LO the two photons 
are constrained to have vanishing total transverse mo- 
mentum Qt = \pn + P2t\, which leads to y* < 0.94. 
At NLO, events with a radiated gluon can have nonzero 
Qt, which removes the constraint on y*. For small Qt, 
the NLO cross section is very unstable and must be re- 
summed in a s In 2 (Qt/pt), as in refs. [2^, Similar 
phenomena have been described in earlier work on iso- 
lated photons |L7|, |j0|. A general description of such 
"edges" has also been given In fig. |l0|, \ti(Qt/pt) 

becomes appreciable only in the two bins centered at 
y* = 0.93 and 0.99. The bins to the left of y* = 0.94 do 
not contain these uncancelled logarithms because the vir- 
tual corrections, with LO kinematics, can contribute and 
cancel them. Of course, all bins to the right of y* = 0.94 
are effectively being calculated at LO, hence their overall 
normalization is not as trustworthy. 

One might be concerned that the parton-level normal- 
ized distributions shown in fig. 10 will be distorted by 
higher-order terms, soft physics, and detector effects, ren- 
dering the information inadequate for improving the sig- 
nificance of the signal. However, this is not the case, 
because 1) only an approximate knowledge of the relative 
shapes of signal and background is required to get most 
of the benefit, and 2) the background distribution can be 
measured experimentally. Once a putative peak is iden- 
tified in the M 77 spectrum, the y* distribution in the 
sideband regions above and below the peak can be mea- 
sured. This information, along with that in fig. [l0[ can 
be used to estimate the true signal distribution, includ- 
ing detector effects, etc. One can then apply an optimal 
observable or maximum likelihood analysis similar to the 
one described above. 

The final distribution that we consider is that of y 77 , 
defined to be the rapidity of p 77 = p 7l + p 72 , the four- 
vector sum of the two photon momenta. For the case 
of the Higgs signal, this is just the rapidity of the Higgs 
boson. We plot the normalized distribution in y 77 in Fig- 
ure [H], for the same choice of mass bin, isolation cuts, 
and scale choices as for Figure |l(]. The difference be- 
tween the signal and the background distributions can be 
mostly attributed to the different parton luminosities in- 
volved in the production; the Higgs signal is produced by 
a predominantly gg initial state, whereas the di-photon 
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FIG. 11: The signal and background distributions as a func- 
tion of 2/ 77 for 116 GeV < M 77 < 120 GeV. The renormaliza- 
tion and factorization scales are fiR — fiF ~ 0.5M 77 and the 
smooth cone isolation cuts are R = 1, e = 1. The Higgs mass 
is taken to be ran = 118 GeV. 



background gets significant contributions from each of qq, 
qg, and gg. In this case the use of y 77 in a maximum like- 
lihood function analysis would improve the significance 
by less than a percent. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we presented a next-to-leading order 
study of the irreducible di-photon background, including 
the corrections to the gluon fusion subprocess gg — > 77. 
The NLO gluon fusion is the largest of the higher order 
contributions not included in previous studies @, 0. 
The scale dependence of the gluon-fusion contribution at 
NLO is roughly the same as at LO if the renormalization 
and factorization scales are varied independently, but is 
significantly reduced if they are varied in unison. More- 
over, the NLO corrections to the gluon-fusion component, 
as a contribution to the total irreducible background, 
were found to be modest, suggesting that this calculation 
is under adequate theoretical control. Indeed, the NLO 
K factor for the gg — > 77 subprocess is only about 65% 
of the NLO K factor for Higgs production. Experimental 
studies using the latter K factor to estimate the former 
one jn| have therefore been a bit too conservative. 

Using the improved calculation we investigated the sta- 
tistical significance of Higgs production as a function of 
the photon isolation cuts. We found that the significance 
depends only weakly on the isolation cuts. Although we 
found a slight enhancement with more severe cuts, we 
noted that as isolation becomes tight, instrumental noise, 
soft hadrons and overlapping events can render the cuts 
experimentally unworkable. Moreover, the perturbative 



predictions become subject to large corrections and are 
unreliable. A better procedure is to include a veto on 
jets near the photon candidates. This can suppress the 
background without suffering from the drawbacks of tight 
photon isolation. As an example, we computed the extra 
suppression due to a jet veto when the effects of frag- 
mentation can be neglected. We found that although the 
background is suppressed, the statistical significance is 
hardly altered. More generally one would need to also 
include a jet veto on the fragmentation contribution, but 
only a weak dependence may be anticipated. The most 
robust improvement we found in the statistical signifi- 
cance of the Higgs signal, albeit still modest, was ob- 
tained using the rapidity difference y* distribution of the 
decay photons. It would be interesting to explore other 
variables characterizing the distribution of hadronic en- 
ergy in the events; a strategy which optimizes the use 
of this information without cutting out large numbers of 
events may be possible. 

These results would need to be investigated further in 
a more realistic simulation than just the parton-level one 
we have done here. In particular, the effects of instrumen- 
tal noise and overlapping events must be included |fl3| . 
One would also need to include a detailed study of the 
reducible tt° background contributions 
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. Once 

a more realistic study is set up, the entire range up to 
rriH < 140 GeV would need to be investigated, instead 
of j ust th e single choice of tuh = 118 GeV used in sec- 
tion IV Q When the LHC di-photon data becomes avail- 
able, the information provided by the sideband regions 
will also be crucial. 

As described in the introduction, there are a number 
of O(o.g) contributions that still have not been included. 
These corrections are all expected to be smaller than the 
O(a^) gluon fusion contribution incorporated into the 
computation presented here. Nevertheless, for complete- 
ness as well as to confirm that there are no surprises, it 
would be useful to evaluate all remaining 0(a 2 s ) contri- 
butions. It would also be useful to incorporate a resum- 
mation of the large logarithms which appear at the kine- 
matic edges of angular distributions |Q and for small 
cone sizes Q Q . 

We are hopeful that further studies will lead to a better 
understanding of the di-photon background, and to an 
increased sensitivity for the Higgs search at the LHC. 
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